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Abstract — In this paper, we focus on two-dimensional connec- 
tivity in sparse vehicular ad hoc networks (VANETs). In this 
respect, we find thresholds for the arrival rates of vehicles at 
entrances of a block of streets such that the connectivity is 
guaranteed for any desired probability. To this end, we exploit 
a mobility model recently proposed for sparse VANETs, based 
on BCMP open queueing networks and solve the related traffic 
equations to find the traffic characteristics of each street and 
use the results to compute the exact probability of connectivity 
along these streets. Then, we use the results from percolation 
theory and the proposed fast algorithms for evaluation of bond 
percolation problem in a random graph corresponding to the 
block of the streets. We then find sufficiently accurate two- 
dimensional connectivity-related parameters, such as the average 
number of intersections connected to each other and the size 
of the largest set of inter-connected intersections. We have also 
proposed lower bounds for the case of heterogeneous network 
with two transmission ranges. In the last part of the paper, we 
apply our method to several numerical examples and confirm 
our results by simulations. 

Index Terms — Connectivity, VANET, queueing networks, spa- 
tial traffic distribution, percolation theory. 

I. Introduction 

Vehicular ad hoc networks (VANETs) play a key role in 
future for safe transportation and easy communications among 
the riders of vehicles. Compared to other types of mobile 
ad hoc networks (MANETs), VANETs have some distinct 
features, i.e., roadmap-restricted mobility patterns as well as 
higher-speed mobility of communications nodes. 

Several papers have considered VANETs at different aspects 
[1]. Among them, connectivity is of crucial importance. The 
main reason is due to strong dependency of network efficiency 
and routing algorithms on this aspect. In fact, in a VANET, we 
have the option of carry and forward approach [2], as well as 
radio propagation for packet transfer. Since radio propagation 
is obviously very faster than carry and forward approach on 
one hand, and the transit delay is of crucial importance in 
achieving the basic goal of VANETs, i.e., safe transportation, 
on the other hand, connectivity plays the main role in routing 
decisions. Succinctly, by knowing the connectivity status at 
different intersections, routing decisions can be made in a 
manner in order to attain minimum packet transfer delay. 

Up to now, many papers have been focused on connectivity 
of MANETs. The authors in [3] have considered the connec- 
tivity of an ad hoc network with a finite number of nodes 



uniformly distributed in a one-dimensional network. An exact 
formula for probability of connectivity for the above network 
(one-dimensional connectivity), can be found in [4]. In [5], the 
authors considered one-dimensional random waypoint mobil- 
ity model for the nodes and then obtained network connec- 
tivity. Some of the other researches have been concentrated 
on the effect of transmission range on connectivity in both 
homogeneous and heterogeneous situations [6]. Moreover, the 
authors in [7] have addressed the effect of random power 
selection against fixed power selection in tradeoff between 
connectivity and interference. In this respect, [8] has also been 
focused on the relation between connectivity and interference 
in a general mobile ad hoc network with spatial Poisson 
process, based on percolation theory. 

Actually, in this paper, two nodes are connected provided 
that they are in the transmission range of each other and 
signal to interference ratio exceeds a threshold. So, the results 
in [8] are very useful when CDMA is the multiple access 
technique in the network. Recently, a performance evaluation 
of connectivity for a VANET has been considered in [9] along 
a highway. And loose lower and upper bounds for connectivity 
along a street have been reported in [4], recently. 

To the best of our knowledge, the papers focusing on 
VANET connectivity including the above mentioned ones, 
only consider one-dimensional connectivity, i.e., along a street 
or highway. But for a typical block of streets including 
several intersections we need two-dimensional connectivity. In 
fact, there are many paths between two typical intersections 
several streets away from each other. When the number of 
streets between two intersections increases, the connection 
possibilities increases more and more. Hence, connectivity 
between two typical intersections becomes more complex. 

Another important problem in VANETs is how the behavior 
of vehicles at different intersections and the arrival rate of 
vehicles at entrances affect the connectivity. Depending on 
the resulting arrival rates of vehicles along different streets 
at both directions, we may have distinct probability of con- 
nectivity along different streets. Proposing a computational 
framework in order to find the relation of these parameters 
and connectivity has not been reported in the literature. In 
this paper, we will fill the above gaps. To this end, we exploit 
a recently proposed mobility model [4] for VANETs. This 
mobility model is based on a BCMP queueing network, i.e.. 
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a quasi-reversible queueing network with nice product-form 
solution property. By solving the related traffic equations we 
are able to find the spatial traffic distribution and especially the 
arrival rate of vehicles at both directions of each street. Then, 
we find the exact probability of connectivity along each street. 
Afterwards, we model the streets and intersections as a lattice 
random graph such that each two consecutive intersections 
has an edge (bond) with a probability corresponding to the 
probability of connectivity of that street. Then, with resort 
to percolation theory and fast computational methods for 
bond percolation proposed in the literature, we are able to 
obtain the related results of connectivity among intersections. 
However, for the case in which we have asymmetric situations 
among different streets, we obtain lower and upper bounds 
for connectivity-related parameters. At all cases, we also 
obtain a threshold for vehicles arrival rate at each entrance 
in order to have the desired connectivity-related results. We 
will see that the general behavior and the results predicted 
by applying percolation theorems on infinite-size lattices are 
also applicable to finite size VANETs. At the last part of the 
paper, we consider several numerical results obtained by our 
approach and confirm our results by simulations. 

Following this introduction, we review the mobility model 
exploited in this paper, in Section HI] and find the spatial 
probability distributions based on vehicle's arrival rate and 
the behavior of vehicles at intersections. In Section |III] we 
compute the probability of connectivity at a typical street with 
respect to spatial distribution at both directions. Afterwards 
we will establish a bound for the case of heterogenous trans- 
mission ranges. Then, we make a brief review on percolation 
results in Section |IV] We discuss fast computational method 
for bond percolation, proposed in the literature in section |V] 
In Section IVII we illustrate the efficiency of our approach 
in a number of scenarios for both symmetric and asymmetric 
situations. We also confirm our numerical results by simulation 
in the same section. We conclude this paper in Section IVIII 

II. Review On Mobility Model 

As discussed in the previous section, we have exploited 
a mobility model recently proposed in the literature [4]. 
Although two mobility models in [4] have been proposed 
for sparse and dense situations, respectively, we only focus 
on sparse situations since in dense situations, there is no 
problem for connectivity. In other words, the dense streets 
are connected with probabilities almost equal to one, so they 
do not require any special concern. Thus, in what follows, 
we discuss the mobility patterns and their modeling, only for 
sparse scenarios 

Obviously, the mobility pattern of each vehicle is strongly 
affected by the topology of the roads. So, the first step for 
considering a mobility pattern, reverts to the road topology 
at the desired region. In this paper, we consider a region 
similar to the one illustrated in Fig. 1, however, there is not 
any restriction on the shape and direction of the streets. In 
the desired topology (Fig. 1), we have several intersections 




Fig. 1 . A typical roadmap and tlie method of mapping different parts of the 
streets onto equivalent queueing nodes. 



and streets. We assume that the vehicles arrive to the desired 
region from entrances according to a Poisson process. Since 
the behaviors of vehicles along the streets are usually different, 
we differentiate the movement of a vehicle along different 
parts of a typical street. In this respect, we make a general 
assumption; a typical vehicle at each street has three phases 
of movement. 

The first phase of movement indicates a transient phase that 
corresponds to front part of the street. In this phase the vehicle 
has a changing speed until it reaches to a nearly constant 
speed. The vehicle maintains its speed at the middle part of the 
street, equivalent to the second phase of its movement. This 
constant speed depends on the driving habits of the driver and 
other factors in the street. By considering the vehicles at that 
street in a statistical sense, we have assumed three categories 
of speeds; low, medium, and fast. Each category at the above 
division has a range of speeds as well as a general distribution 
(e.g., uniform), such that each vehicle of a typical category 
selects a speed from the corresponding range according to the 
corresponding distribution. After moving along the middle part 
of the street, we have considered another phase of movement 
due to approaching to an intersection. We have considered 
two categories of speeds in general for this end part of the 
streets; low, and medium. Low speeds represent the situation 
of encountering stoplights, stop-signs, or any other means of 
blocking. And medium speed category corresponds to the case 
of not encountering an effective obstacle. Correspondingly, we 
have considered two speed categories for the vehicles at the 
front part; one category represents the vehicles arriving at the 
front part with low speed and the other category represents 
the vehicles arriving at the front part with medium speed. 
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Therefore, we have two speed categories at both front and 
end parts of each street and three speed categories for the 
vehicles at the middle part of the streets. Usually, we assume 
that the length of the front and end parts of the streets are 
very smaller than the length of the middle part. It is worth 
mentioning that each speed category represents a range of 
speeds with a specific general distributions (e.g., uniform, 
truncated normal, etc.). 

For the sake of simplicity, in this paper we have assumed 
that through a street a vehicle does not change its direction 
except at the intersections. Also, a vehicle at an intersection 
can go straight forward, turn to the right or turn to the 
left. However, these assumptions do not Umit the considered 
mobiUty pattern. 

In the above discussed mobiUty pattern, if the density of the 
vehicles at a typical street is low, i.e., a sparse scenario, we can 
consider the mobihty of the vehicles independent of each other. 
Actually, in this scenario when a vehicle approaches to the 
slower vehicle at its front, the faster one overtakes the slower 
one with a simple maneuver. Therefore, the independency 
assumption is correct providing that such a maneuver is 
possible. Clearly, this is the usual case for sparse scenarios. 

Now, we consider a queueing network comprising several 
nodes and a number of customers with several classes. In our 
modeling approach, we map the vehicles onto the customers 
of the queueing network. And we consider three nodes for 
each street at each direction, corresponding to front part, 
end part, and middle part of that street (Fig. 1). Also, we 
correspond the residence time of a vehicle at each part of the 
street to the service time of the corresponding customer at 
the corresponding node. Therefore, in the queueing network, 
when a customer departs from a node representing the front 
or middle part of a typical street it is routed definitely to a 
specific node, i.e., the node representing the middle or end part 
of the same street, respectively. However, when a customer 
departs the node representing the end part of a street, i.e., the 
corresponding vehicle reaches to an intersection, it is routed to 
three nodes representing the front parts of other three streets 
at that intersection with some probabiUties (Fig. 1). 

As we discussed we consider different speed categories at 
different parts of the streets, representing approximately dif- 
ferent behavior of the vehicles along a street. In the proposed 
mobility model we simply map each speed category onto a 
customer class. Therefore, we have three customer classes at 
nodes representing the middle parts of the streets, and two 
customer classes for the nodes representing the front and end 
parts of the streets. Each customer class at each node has a 
random service time with an arbitrary distribution equivalent 
to the distribution of the corresponding speed category at 
the corresponding street. Then, the average service time of 
a customer class at each node is determined with respect to 
the distribution of its corresponding speed. 

Therefore, we consider each node at the above queueing 
network, as an M/G/oo node. So, in the case all parts of 
the street are in sparse situation, we actually have a queueing 
network comprised of M/G/oo nodes with 2 or 3 customer 



classes, corresponding to middle or front/end parts, respec- 
tively. This is a BCMP queueing network with a product-form 
solution. In this network the following traffic equations are 
satisfied [10]: 

N C(k) 

aju = ^ OikvrkvjuJ = 'i-,---,N;u=l,...,C{j) 

(1) 

where Xju is the exogenous arrival rate of class-u customers 
at node j, C{j) is the number of classes at node j, and 
fkvju denotes the routing probability of a departed class- 
V customer from node k to node j as a class-u customer. 
Moreover, Aj„ equals zero for the internal streets (see Fig. 1). 
Also, Xju in (1) denotes the arrival rate of class-w customers 
at node j in the network, representing exogenous arrival rates 
as well as routing from other nodes. It is obviously equals the 
departure rate of class-u customers from node j in the case 
of stability. We also have the following relation for routing 
probabilities [10]: 

'^'^rju,kv+rjufi = l;j = l,...,N;u= (2) 

where denotes the world outside the network. Then, we 
obtain the spatial traffic distribution corresponding to the 
proposed BCMP queueing network comprised of M/G/oo 
nodes, as in the following [10]: 

N 

7r(n) = ]j7r^-(n^-); 

n = (ni, n2, . . . , un), nj = (n^i, . . . , rijcu)), (3) 

CU) n,^ ^ C{j) 

Where denotes the inverse of average residence time, 
i.e., service rate, of a class-M customer at node j, and C(j) 
denotes the number of classes at node j. In (3), 7r(n) represents 
the probability of nn class-1 customers at node 1, ni2 class- 
2 customers at node 1, etc. By (3) we obtain the probabiUty 
of existing different number of customers of different speed 
categories at different parts of the streets, i.e., spatial traffic 
distribution. Then, by solving the traffic equations and finding 
the arrival rates of vehicles at each street, we are able to find 
the probability corresponding to different number of vehicles 
at that street. We employ this result in the next section for 
computing the probability of connectivity. 

III. Analysis of Connectivity 

In this section we compute the probability of connectivity 
of a street. As it is mentioned before our city is modeled by a 
square graph and each street is divided to three sections. And 
also, it turns out from previous discussions that the number 
of nodes at each section is a Poisson random variable. Let 
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Parameter 



Low speed in the middle part 
Medium speed in the middle part 
Fast speed in the middle part 
Low speed in the front part 
Medivmi speed in the front part 

Low speed in the end part 
Medivmi speed in the end part 
Length of the middle part 
Length of the front part 
Length of the end part 
Transmission range 



Value 



Uniform [3, 14] rajs 
Uniform [14, 22] mjs 
Uniform [22, 33] mjs 
Uniform [0.3, 3] m/s 
Uniform [3, 14] mjs 
Uniform [0.3, 1.5] m/s 
Uniform [1.5, 14] m/s 

1600 m 

200 m 

200 m 

200 m 



TABLE I 

Typical values used for parameters of our analysis 
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Fig. 2. Four possible states for inserting r nodes of type 2 into q blocks. 



pi, p2 and p3 denote the parameters of the Poisson random 
variables corresponding to a typical three sections of the street. 
Since the length of first and last sections are comparable with 
transmission ranges (see table I), it is just necessary to have 
at least one node in these sections to be sure of connectivity. 
We state similar arguments for connectivity of middle section 
as in [11] and we extend it, in order to find a lower bound 
for connectivity corresponding to heterogenous transmission 
ranges. Thus, we have the following relations: 

P{con) = P{con_seci) x P{con_sec2) x P{con_sec3) 
= Pseci {Nodes 7^ 0) X P{con_sec2) x PsecsiNodes ^ 0) 
(1 - e-P') X P{con_sec2) x (1 - e"''^), 

(4) 

where P{con), P{con_seci), Pgec, (Nodes ^ 0) denote the 
probabihty of connectivity along the street, probabihty of 
connectivity at section i of the street and probability that 
section i contains at least one vehicle, respectively. Using 
the results of [11] for one-dimensional connectivity in and 
considering that the length of the second section of the typical 
street and transmission ranges are D and R, respectively, we 
conclude that: 



N 



P{con_sec2) = ^ P{Nodes = n) x P{con\Nodes=n) 



n=0 



n=0 i=0 ^ ' 



(5) 



where m = min{n+ 1, [-^J}, for each n in the summation. 

Let us assume that we have a street of length D and every 
vehicle has a transmission range equal to .ti and X2 (.xi < X2), 
with the probabihties p and (1 — p), respectively. We mention 
nodes having x\ as a transmission range as nodes of "typ^ 
1" and others as nodes of "type 2". Besides, let (") be zero 
wherever i > n or z < in the following discussions. 

Theorem 3.1: Suppose we have N nodes distributed in- 
dependently and uniformly in the street. The probabihty of 
connectivity in the street is lower bounded by: 



^p^{l-p)''-^Q{r,N), 



(6) 



r=0 



where, 

Q{r,N) = 



N 



(7) 



for r > 0, in addition Q{0, N) = p(0, -h 1) and 

p(ni,n2) = 



Af+l 



where the summation is taken over positive values of 1 — Ixi — 

(k - l)X2. 

Proof: Suppose that Ui represent the distance of i*'* node 
from the origin. Ordering these random variables (RVs) we 
obtain: 

< [/(I) < . . . < f/(iv) < D, (9) 

Define: 

^ TT TT TT ^ n TT ^ 



= ?7(i) - {7(i_l), t/'(0) = 0, ?7(JV+1) = D, 

Wi = l,2,...,N+l. 



(10) 



Again, considering the results for street connectivity [11], 
for c,^ > 0, {a^ IJ^o^C {1, 2, . . . , iV} and E^^i c., < ^ we 
know that (1— X^^^i -^)^ expresses the probability -PjAjj > 
Cti,...,Ai^ > Ci^.}. Now, using principle of inclusion and 
exclusion we attain the following relation: 



r 

= E(-i)-'' E p{\.>ci„...,\^>c,j 



3=0 



(11) 



Ei-^y E (i-Et)" 



j=0 1=1 



We take this formula as our template for the corresponding 
calculation. We shall distinguish nodes with different ranges 
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and their positions in order to find our convenient result. For a 
moment, we restrict ourselves for a moment to a specific case 
in which nodes having X2 as transmission range are arranged 
in a way that there are only q blocks of consequent nodes of 
type 2. And also assume that only r nodes are of type 2. It is 
worth mentioning that this case happens with the probability 
p'^{l — p)^~^. By previous assumption, there are only four 
cases for locating these q blocks and Fig. |2] depicts them. In 
case 1 and 2, nodes of type 1 are also located in q blocks 
whereas they are located in g + 1 and q — 1 blocks in respect 
to case 3 and 4. 

In Fig. |2] Xi denotes a block of consequent nodes of 
transmission range Xi. Notice that each state has just q blocks 
including nodes of type 2. In order to compute the number of 
situations related to case 1 in Fig. |2] we have to put r nodes 
of type 2 into q blocks and other N — r nodes of type 1 into 
other q blocks which is given by the equation (^^^) (^I^Ig)- 

To be sure of connectivity, it suffices for the distance of 
nearby nodes (A^s) to be less than the transmission ranges 
of those nodes. And also, in order to have connectivity with 
the entrance and exit part of the street, Ai and Aat+i should 
be less than and X(^j^y respectively. But there also exists 
cases in which street is connected and the distance between 
some consequent nodes is greater than one of the nearby 
nodes' transmission range. For example, this can occur when 
a node of type 1, lying between two connected nodes of type 
2, is isolated from them. Thus, as we do not consider these 
cases, our assertion is a lower bound for connectivity. Next, 
we choose c^s in (fTTI) to be a; i or X2 so that our aim will be 
achieved. Assertion (fTTT i is symmetric and linear, thus it only 
depends on the number of c^s which are equal to Xi. In case 
1 , among A^ s (r + q) distances are X2 and other N+l — r — q 
ones are xi. This symmetric property and similar arguments 
for other cases shed light into the following equation (after 
replacing by xi or X2, (fTTT i will change to (O): 



^^•^(1 -p)^-''Q(r,7V) < P(con), (12) 



where, 

Qir,N) 



N 



N 



(13) 

for r > and Q{0, N) = p{0, N + 1). One can ignore the 
differences between ^(A^+l—r+g, r—q),p{N—r+q, r—q+1) 
and p{N ~ 1 ~ r + q,r — q + 2) and approximate (fT3] l by 
dTji. This approximation is on the behalf of obtaining a lower 
bound. And this completes our proof. 
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Fig. 3. Probability of bond percolation 



It is clear that the inequality reduces to equality for p = 
0, i.e., one-dimensional connectivity with fixed transmission 
range. 

IV. A Review On Percolation Theory 

Percolation theory was originally developed by Hammersley 
and Broadbent [12] in 1957 as a model for describing fluid 
spread through porous mediums. It basically deals with the for- 
mation of long-range connectivity in random systems. Later, 
it found many applications in a wide range of phenomena 
that exhibit phase transitions and critical behavior especially 
in statistical physics and materials science. More recently, 
concepts of percolation theory have been applied in the study 
of the network topology, connectivity and robustness to yield 
valuable results. 

Consider an infinite 2D square lattice grid in which each 
of vertices or edges are independently occupied (open) with 
probability p and empty (close) with probability 1 — p. The 
former case is called site percolation and the latter is known 
as bond percolation. We will restrict our discussion to bond 
percolation problem, however, similar arguments hold for the 
site percolation case as well. 

For the bond percolation case two open bonds are con- 
nected if they share a common vertex or there exists a 
path of subsequent open bonds between them. Clusters are 
then formed through connections between nearest neighbor 
occupied bonds. The network is said to percolate if the size 
of cluster containing a given vertex, e.g., origin, is infinite. 
More formally the percolation probability function is defined 
as follows [13,14]: 

0{p) - Prp{\C\ = ^) 

where C is the cluster containing the origin. 

Fig. 3 shows a typical plot of the percolation probability. It 
is a fundamental result of percolation theory that there exists 
a critical probability threshold [13] for which 

f 0{p) = for p < 

\ 9(jp) > for p> pc 

where, 

< = sup{p : 6'p = 0} < 1 
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The value of critical probability depends on the lattice struc- 
ture, dimensionality of the grid and type of percolation system. 
Unfortunately the exact values of the percolation thresholds are 
known for only a few special cases. The following non-trivial 
theorem for bond percolation threshold in square lattices wiU 
be enough for our purpose. 

Theorem 4.1: The critical probability of bond percolation 
on the square lattice is 1/2. 

Another important result from percolation theory is the zero- 
one law for the existence of infinite open cluster and its 
uniqueness. 

Theorem 4.2: The Probability that an infinite open cluster 
exists is if p < Pc and 1 if p > Pc- The infinite cluster is 
almost surely unique whenever it exists. 

This infinite cluster is also called spanning cluster since it 
spans the entire lattice. When the value of bond probability 
is relatively small, the grid mainly consists of isolated islands 
of nodes with a few pairs of neighboring bonds. As the bond 
probability increases, these islands grow and some of them 
merge to form larger clusters. When the occupancy probability 
is about to exceed the critical threshold, suddenly all clusters 
are absorbed into one dominating cluster and the unique 
infinite cluster appears. Above this threshold, other small size 
clusters gradually join the spanning cluster. So clearly is 
exactly zero in the subcritical phase, since there is no spanning 
cluster in this phase. 

For the finite size lattices, however, instead of a single 
point there exists a transition region where the percolation 
probability raises sharply with the increase in the value of 
edge probabihty. Here the transition region becomes narrower 
with the growth of lattice size. We wiU discuss this matter 
later in numerical results. 

V. Computational Algorithm 

In this section we will briefly introduce an efficient Monte- 
Carlo approach estimating desired parameters in percolation 
systems. It is clear that testing all permutations of a bond 
percolation problem in a N x N grid is computationally 
intractable, since there exists 2N{N — 1) bonds and that the 
solution requires 0(22^(^-1)) of time even for a single value 
of edge probability. Therefore, using a Monte Carlo approach 
seems to be inevitable. 

Newman and Ziff [15] proposed a fast algorithm for the 
estimation of percolation probabilities and lattice quantities 
directly for any arbitrary value of edge probability. The idea 
is that in order to estimate any parameter, e.g. average or 
giant cluster size, in a finite size percolation it suffices to 
only estimate that parameter for all possible numbers of total 
edges in the corresponding lattice and then use the law of total 
probability: 

Q{p) =Y,P{M = m)Qm = J2 

m m ^ ^ 

(14) 



where M denotes the total number of edges in the lattice and 
Qm represents the expected value of the desired parameter in 
a lattice with m total occupied bonds. 

Now the problem reduces to the computation of Qm coeffi- 
cients. Consider for example "average cluster size" and "giant 
cluster size" as two parameters of interest. So we need to find 
contiguous clusters of open bonds. One simple suggestion for 
the estimation of these parameters is to independently generate 
a number of configurations of the percolation system over all 
possible states with total M edges. The clusters can then be 
obtained using a simple depth-first search (DFS) or breadth- 
first Search (BFS). However one can easily check that this 
approach is not computationally effective and is intractable 
for large values of lattice size [15]. 

We can do much better by noting that any state with n + 
1 bonds is generated by adding an extra randomly chosen 
bond to a sample state with n bonds. Thus we can start with 
an empty lattice and add one bond to it at each step. The 
advantage is that we can keep track of the clusters and update 
them at each step within much less operations since each bond 
can merge two clusters at most. 

By repeating this procedure for a number of iterations and 
averaging over the sample lattice quantities we expect to have 
a good estimate of these values. This is because the precision 
of a statistical measurement depends on the number of samples 
taken and not the size of the sample space. 

To keep track of the clusters and update them at each step 
we need to find an efficient data structure to perform two 
operations efficiently. These operations are Find(x) that returns 
the cluster that node x belongs and Union (x, y) to merge 
two clusters that contain x, y nodes whenever we add a bond 
between them. This problem has been extensively studied in 
computer science and there exists a family of efficient algo- 
rithms for it, generally known as Union-Find Algorithms. We 
can summarize this algorithm for bond percolation problem in 
an N X N square grid in the following steps: 

• Initialization: Consider each of N'^ nodes in an empty 
lattice as single clusters. Each node will have a pointer 
to its parents, initially all the pointers will be null. 

• Find(x): We start from x, and recursively traverse up the 
tree following the parent pointers till we get to the root 
which will be a node with null parent pointer. Return 
the root as the cluster which contains x. Each time this 
procedure is called we set all the parent pointers of all 
the nodes on the path to the root directly to the root of 
the tree (Path compression step). 

• Union(x, y): First we perform two find operations to 
obtain the corresponding cluster roots of each node. If 
both nodes belong to a single cluster we do nothing, 
otherwise we merge them by setting the parent pointer 
of smaller cluster root to the parent pointer of the larger 
cluster's root. Thus we should anly keep track of the 
cluster sizes for each of the roots. 

Analysis of this algorithm [15] shows that performing a 
sequence of m operations over n elements requires an overall 
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Fig. 4. Connectivity properties of finite bond percolation as a function of edge probability, for different square sizes. 
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Fig. 5. Connectivity properties of VANETs as a function of entrance rate of vehicles, for different city sizes. 



running time of 0{{n+'m)log* (n)) which grows almost linear 
with the size of grid nodes in practice. It should be noted 
that this algorithm can be easily modified to find percolation 
thresholds and different properties for other lattice shapes and 
percolation systems as well. 

VI. Numerical Results and simulations 

In this section, we will investigate our approach to analyse 
connectivity features of the VANETs in a number of scenarios. 
We will consider three measures for the connectivity: "Giant 
Cluster Size" which is the fraction of the intersections that 
belong to the largest connected cluster, "Average size of 
connected clusters" or simply "Average Cluster Size" and fi- 
nally "Probability of Perfect Connectivity", i.e., the probability 
that all intersections are connected to each other In all our 
scenarios, we will use typical values shown in Table I, for 
different parameters of the presented mobility model (section 

Scenario 1. Connectivity as a function of Edge Probability 

First, we examine connectivity properties of the finite bond 
percolation problem as a function of edge probability (p). In 
our model of the city, this will correspond to the connectivity 
of the intersections of the city as a function of the street 
connectivity. These parameters are illustrated in Fig. |4] for 
different lattice sizes. The emergence of giant cluster size 



at about the edge probability of 1/2 completely agrees with 
theorem [4. 1 1 (Fig pl-(a)| i. As expected, this transition is sharper 
for larger values of lattice size. In Fig. |4(b)[ |4(c)[ we observe 
that the transition for average cluster size and probability of 
perfect connectivity occurs in a larger value of edge probability 
depending on the lattice size. This is due to the fact that, even 
after the appearance of a unique giant cluster (theorem 14.21 ). 
there are still many isolated islands of nodes in the lattice. 
It is not also surprising to see the increase of the transition 
probability for larger values of lattice size. 

Scenario 2. Connectivity as a function of Entrance Rate 

We demonstrate the results of our analysis of connectivity 
properties as a function of the entrance rate of the vehicles. 
We use our queueing network model to find the probability of 
street connectivity and then apply the fast algorithm proposed 
in section [V] to find the desired parameters. Fig. |6] shows the 
probability of a street connectivity as a function of entrance 
rate of the vehicles (|4]i. The connectivity properties, 
previously mentioned, are demonstrated in Fig. |5] 

It is also necessary to validate the analytical results of our 
method with the aid of a computer simulation. With the same 
parameter setup, the simulation results are compared with 
analysis results in Fig. [T] The congruence of these diagrams 
is approved our method of connectivity analytical for 2D 
vehicular networks. 
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Fig. 8. Connectivity propeilies of VANETs as a function of transmission range of Veliicles, for different city sizes. 



Scenario 3. Effect of Transmission Range 

As an application of our method, we now consider the effect 
of the transmission range of the vehicles on the connectivity 
properties of the network. The problem of finding optimal 
range for radio transmission is one of the main concerns of 
the wireless design engineers, since choosing an inappropriate 
transmission range, whether low or high, is destructive for 
the network. The disadvantage of low-power transmission is 
that we obviously lose the network connectivity. On the other 
hand, a very high power of transmission is not only a waste 
of energy, but also a cause of interference among nodes in 
the network. Fig. [8] displays the connectivity properties of the 
VANETs with respect to the transmission ranges of the vehi- 
cles for several values of arrival ranges, i.e., A — 0.1,0.2,0.3 
respectively. Clearly, given an estimation of entrance rate of a 
city, we can find the optimal value for the transmission range 
of vehicles. 

Scenario 4. VANET Connectivity in an asymmetric city 

In our mobility model, we considered a symmetric model 
for the routing of vehicles between distinct intersections. 
However, in a real situation this assumption is not valid 
because different parts of the city have certainly different 
loads of traffic. We can simply generalize our mobility model 
to solve this problem by assigning a traffic weight factor 
to each intersection of the city. The vehicles will then be 
routed between intersections according to this factor, i.e., the 
probability that a vahicle chooses one intersection as its next 



destination will be proportional to the traffic weight of the 
corresponding intersection. 

In this case, the street connectivity probabilities of the city 
are not equal any more. However, note that we can find a lower 
and upper bound for connectivity characteristics by assuming 
that the probability of all streets equal to the least and largest 
probabilities among all streets of the city, respectively. As 
an example, we consider a city in which all traffic weights 
are distributed uniformly in [1,2]. Fig. |9] illustrates these 
bounds for a sample city with such a traffic distribution model 
mentioned above. We are much more interested in the lower 
bounds, since upper entrance rates will guarantee the network 
connectivity. The upper bounds gives us a rough estimation of 
the precision of our lower bound. 

Scenario 5. Multiple Transmission Ranges 

Finally, we take our theoretical results of section |lll] into 
account for the problem of street connectivity in a network 
of vehicles with two transmission ranges. Here, each node is 
given transmission range 200m or 400r7i with probability p 
or (1 — p) respectively. Obviously, every street connected by 
nodes having transmission range less than or equal to 2R will 
be connected if we set all transmission ranges to 2R. Thus, the 
graph related to p — can be thought of as a upper bound for 
other graphs. Fig. [TO] shows lower and upper bounds provided 
by Theorem 13. II 
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Fig. 9. Bounds on connectivity properties obtained from minimum and maximum probabilities of street connectivity in a random asymmetric city. 
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Fig. 10. Lower bounds on connectivity properties (section Hill for different values of p, the probability that a node takes xi = i? as its transmission range. 
(X2 = 2R,R = 200m). 



VII. Conclusions 

In this paper, we considered two dimensional connectiv- 
ity of VANETs and proposed a new method for evaluating 
connectivity measures such as average size of the connected 
intersections and size of the giant cluster of them, as a function 
of the entrance rates of vehicles. Our approach was based 
on modeling the mobility of the vehicles with a BCMP open 
queueing network to find spatial distribution of the vehicles 
and exact probability of connectivity for each street. We then 
exploit results of percolation theory and a fast algorithm 
for computing and interpreting connectivity properties of the 
large scale network. We also derived lower bounds on the 
connectivity probability of each street for the heterogeneous 
case in which vehicles could have two different transmission 
ranges. Finally we illustrated the efficiency of our approach in 
a number of scenarios, each yielding to some valuable results. 
These included finding bounds on connectivity features in a 
city with asymmetric traffic load and determining optimal 
value for the transmission range of vehicles for a given 
entrance rate of a city. 
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